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By George E. Castile and Robert W. Herr 

SUMMARY 


Results of sixty-six experimental flutter tests are presented to 
show some effects of density and Mach number upon flutter speed. Two 
uniform, cantilever wing models, one having a small amount of sweep, 
were fluttered through a range of Mach number up to 0.82. The mass 
ratio, k, for these models ranged from 0.00375 to 0.11, while Reynolds 

number varied from 0.2 66 x 10^ to. 3.000 x 10^. Two testing mediums, 
having different velocities of sound, were used to obtain flutter at 
two Mach numbers at the same density. 

Multiplying the flutter speeds by a simple compressibility factor 
superimposed the plots of flutter speeds obtained in the two gases 
throughout the range of densities, and this correction was indicated to 
be adequate for the range of parameters considered. An additional 
correction to the flutter speed was made to account for the change due 
to the reduced-wave -length parameter. Correcting the experimental 
flutter speed for both Mach number and wave-length effects .resulted in 
a nearly constant dynamic pressure at flutter. No Reynolds number 
effect upon flutter speed was indicated except at values below 

0.400 x 10^. A specific change in tunnel -wall interference had 
no appreciable effect on the flutter speed. 


INTRODUCTION 


Present-day trends toward flight at high speeds and high altitudes 
have led to increasing interest in the flutter speed at low densities 
and high Mach numbers. An experimental investigation has consequently 
been made to determine some of the effects of these parameters on flutter 
speed. These experiments were carried out on two uniform cantilever 
wings, one unswept and one with 15° sweepback. A study of density 
effects was made possible by fluttering the models through a wide range 
of pressure and the Mach number effects were isolated by the use of two 
testing mediums having different velocities of sound. 
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In the course of this investigation, correlation of the experi- 
mental data required consideration of the effects of the reduced-wave- 
length parameter and Reynolds number. In addition, the effect of the 
tunnel wall parallel to the span was studied by changing the distance 
from wall to model. 
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SYMBOLS 

length of model measured along midchord, feet 

half -chord of model normal to leading edge, feet 

nondimensional coordinate axis of rotation relative to 
midchord 

nondimensional. radius of gyration relative to axis of 
rotation 

nondimensional center of gravity relative to axis of 
rotation 

natural frequency of torsional vibrations, cycles per second 

natural frequency of bending vibrations, cycles per second 

flutter frequency, cycles per second 

density of testing medium, slugs per cubic foot 

mass of wing per foot of length 

Mach number 

mass ratio (*£)■ 

flutter velocity, feet per second 

angular flutter frequency, radians per second (2itff.) 
in-phase component of circulation functions (reference 2) 
out-of -phase component of circulation functions (reference 2) 
sweepback angle of midchord line, degrees 
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Subscripts: 

e experimental value 

1 first mode 

2 second mode 


APPARATUS 


These experiments were performed' in the Langley flutter tunnel. 

The tunnel is of the closed-throat, single -return type in which either 
Freon-12, air, or a mixture of the two may be used as a testing medium. 
A characteristic of this tunnel is that Mach number and density of the 
testing medium can be varied independently. The pressure can be varied 
from atmospheric to 4 inches of mercury absolute, and a density range 
from O.OOO317 slug per cubic foot in air to 0.0106 slug per cubic 
foot in Freon-12 is thus provided. The Mach number range of the tunnel 
is from 0 to approximately 0.92. Flutter was attained at Mach numbers 
up to 0.82. The mass ratio, k , for these models ranged from O.OO375 

to 0.11. The Reynolds number range was from 0.2 66 x 10^ 

to 3.000 x 10 6 . 

In order to investigate some of the effects of tunnel -wall 
interference upon flutter speed, two tunnel test sections were employed: 
one of circular cross section 4.5 feet in diameter; the other, the 
^•5 foot section with side walls installed to give a section approxi- 
-mately 2 feet by 4.5 feet as shown in figure 1. 

Data were obtained for two different model wings (figs. 2 and 3). 
The models were built up of solid balsa wood with dural inserts and 
were mounted rigidly at the top of the test section as cantilever beams 
at zero angle of attack. Model A had 15° of sweepback whereas model B 
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had none. Both models incorporated the NACA 16-010 airfoil section. 
The physical properties of these models were as follows: 


Model 


A B 


l, feet 

b, feet 

m, slugs per foot . . 

O 

r a , nondimens ional . 
a, nondimens ional . . 
Xq,, nondimensional 
fhpj cycles per second 

fh2> cycles per second 
cycles per second 
k f nondimensional 


2.167 

4.0 

0.167 

0.333 

0.0099 

0.0272 

0.107 

0.324 

-0.16 

-0.36 

-0.024 

-0.238 

9.4 

4.5 

56 

27 

97 

47 

8.8lp 

12 . 84p 


Two sets of strain gages were attached to each model near the root 
and were connected to a recording oscillograph in order to determine 
the. natural bending and torsional frequencies as well as the flutter 
frequency. 


TEST PROCEDURE 


Flutter occurs very suddenly and is usually destructive unless 
great care is exercised. Since this series of experiments required 
the physical properties of the models to remain unchanged, preservation 
of the models was an important consideration. The airspeed in the 
tunnel was, therefore, brought up gradually until flutter was observed. 
The airspeed was then reduced immediately in order to prevent destruc- 
tion of the model. At the point of flutter, an oscillograph record 
was taken and the tunnel pressure, dynamic pressure, and tunnel temper - 
ature were recorded. The percentage of Freon-12 in the testing medium 
was obtained after each test. The natural bending and torsional 
frequencies of the model were recorded on the oscillograph after each 
test in order to ascertain whether the model had been damaged or 
weakened by flutter. These frequencies did not change a measurable 
amount throughout' the series of tests. 


DISCUSSION OF RESULTS 


Results of sixty-six experimental flutter tests are presented in 
tables I and II. As stated in the introduction, effects of density, 
Mach number, wave length, Reynolds number, and tunnel -wall interference 
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were partly isolated and the effects of these parameters are shown 
in figures 4 to 11 and briefly discussed in the following paragraphs. 
Although density is the principal parameter, corrections for Mach number 
and wave length must be applied to the flutter speeds before a compre- 
hensive picture of the variation of flutter speed with density can be 
gained. 

Mach number .- Many investigators have studied, theoretically, means 
of correcting the air forces acting on an oscillating wing for the 
effects of compressibility. The calculations involved are generally 
quite tedious. Theodor sen and Garrick have suggested the simple 

correction 4jr~- m 2 •as a first approximation (reference 1). The 
re suit s presented in figures 4 and 5 show that the experimental flutter 
speeds are lower in Freon-12 than in air at the same mass ratio (k). 

This result represents, for the most part, a Mach number effect since 
the speed of sound in Freon-12 is about half of the speed of sound in 
air, and the Mach number is therefore higher in the former medium. 

Inasmuch as model A had 15° sweepback, the Mach number component 
perpendicular to the leading edge (M cos A) was used and the correction 

thus becomes \^/l- (M cos A). In figure 6 the compressibility- 
correction for 0° and 15° sweep is plotted as a function of Mach number. 

Figures 7 and- 8 show the data corrected for compressibility. The 
adequacy of the correction is shown by the fact that the flutter speeds 
in the two mediums were superimposed throughout most of the density 
range. 

Reduced-wave -length parameter . - Another factor which affects the 
air forces, and thereby the flutter velocity, results from the 
vibrating wing itself. The circulation, and therefore the lift of an 
oscillating airfoil is known to depend on the reduced-wave -length 
parameter v/bo>. As v/bo> is increased, the air forces on the wing 
are increased and finally approach the steady-state value. Theodorsen 
has expressed this change in lift by means of the so-called F and 
G functions (reference 2). The total effect on the air forces is pro- 
portional to the absolute magnitude of the vector VF^ + G^. For a 
given density, these air forces are proportional to the velocity 
squared. The correction to be applied to the flutter speed is then 

+ G 2 . This factor is plotted as a function of the reduced-wave- 
length parameter in figure 9» An attempt is made in figures 10 and 11 
to extract the effect of the reduced-wave -length parameter by multi- 
plying the flutter speed by the factor \J}/F 2 + q2. 

Density . - Correcting the experimental flutter speeds for both Mach 
number and wave-length effects resulted in a straight-line plot, which 
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if extended) would pass through the origin. This result would seem to 
indicate that, in order to obtain simple bending-torsion flutter from a 
given model, a constant dynamic pressure is required. It should be 
pointed out, however, that this relationship may not be true for mass 
ratios and frequency ratios beyond the range encountered in this 
investigation. 

Reynolds number .- If the experimental data are reduced to the form 
presented in figures 10 and 11, some insight may be gained as to the 
effects of Reynolds number. In figure 10 (model A), at any given 
density ratio, the data obtained in Freon-12 and those obtained in air 
have been nearly superimposed, even though the Reynolds numbers 
differed considerably. The Reynolds number is therefore concluded to 
have little or no effect upon the flutter speed of this model. Three 
flutter points of model B (fig. 11), however, fell considerably below 
the curve predicted by- theory. These three points were made at test 
Reynolds numbers below 0.400 x 10^. Inasmuch as the flutter speed is 
a function of the slope of the lift curve for the oscillating case, 
the Reynolds number might be expected to affect the flutter speed in 
the low Reynolds number range. The dash curve in figure 11 is the 
theoretical speed as determined by the two-dimensional potential-flow 
theory. These values are for incompressible speeds and have been 

multiplied by the function 'Vf 2 + G 2 in order to permit correlation 
with the experimental data. In the low-density range, fbr high values 
of v1/k, the theoretical curve is noted to continue essentially as a 
straight line. The deviation of the experimental data from theory may 
therefore possibly be an effect of Reynolds number. 

Tunnel-wall interference .- In order to interpret any abnormalities 
in the data that might arise from tunnel side-wall interference, two 
test sections of different configuration were employed. One section 
was 4.5 feet in diameter) whereas the other was the same section with 
side walls installed as shown in figure 1. Differences in flutter 
speed under corresponding conditions in the two test sections can be 
seen in figures 4 and 5 » As these differences were of the same order 
of magnitude as the experimental error, the effects of c han ging the 
tunnel configuration were considered negligible. 


CONCLUSIONS 


An experimental investigation conducted of two uniform wings in 
the Langley flutter tunnel to determine the effects of density and Mach 
number on flutter speed indicated the following conclusions: 
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1. A simple compressibility correction was adequate to account for 
Mach number effects. 

2. After the experimental data were corrected for compressibility 
and wave-length effects, the dynamic pressure at flutter was nearly 
constant over a wide range of densities. 

3. No effects attributable to Reynolds number were found through- 
out most of the range tested. The results appear to indicate, however, 
that a Reynolds number effect may occur below a Reynolds number 

of 0.400 x 10 6 . 

4. A spelcific change in tunnel side-wall interference had no 
appreciable effect upon the flutter speed. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., September 27, 1949 
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Model mounted rigidly at top 



Figure 1»- Cross-sectional view of 2- by ij.. 5-foot 

test section. 
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Figure Ij..- Variation of experimental flutter speed with 
the square root of the mass ratio for model A. 
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Figure 7*" Variation of incompressible flutter speed with 
the square root of the mass ratio for model A. 





Figure 8.- Variation of incompressible flutter speed with 
the square root of the mass ratio for model B. 







Figure 9.- Variation of Wf^ + with reduced wave 
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